The nature of the magnetic and structural phase transitions in BaFe2As2 
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We present the results of an investigation of both the magnetic and structural phase transitions 
in a high quality single crystalline sample of the undoped, iron pnictide compound BaFe2As2. Both 
phase transitions are characterized via neutron diffraction measurements which reveal simultaneous, 
continuous magnetic and structural orderings with no evidence of hysteresis, consistent with a single 
second order phase transition. The onset of long-range antiferromagnetic order can be described 



by a simple power law dependence (p{T) 



[l - with /3 = 0.103 ± 0.018; a value near 

the P — 0.125 expected for a two-dimensional Ising system. Biquadratic coupling between the 
structural and magnetic order parameters is also inferred along with evidence of three-dimensional 
critical scattering in this system. 

PACS numbers: 74.70.Dd; 75.25.+z; 75.50.Eo; 75.40.Cx 



I. INTRODUCTION 

The recent discovery of high temperature superconduc- 
tivity (high-Tc) within the iron pnictidcs^ has resulted 
in a torrent of research into the electronic, magnetic, 
and structural behaviours of this new class of super- 
conductors and their isostructural variants. The pnic- 
tide superconductors now constitute a long sought com- 
parator for testing the physics fundamental to the phe- 
nomenon of high Tc superconductivity through compar- 
ison to the widely studied class of high Tc lamellar cop- 
per oxides. The great deal of the interest surrounding 
this new iron pnictide class of superconductors originates 
from their striking parallels and notable differences rela- 
tive to the high Tc cuprates. Chief among the similarities 
are the layered building blocks comprising the fundamen- 
tal structure of these materials and the proximity of an 
antiferromagnetic ordered phase to the superconducting 
state within their respective phase diagrams^'^'^. How- 
ever, it was also quickly uncovered that the parent com- 
pounds of the iron pnictides are poor metals, and studies 
of magnetism in these materials seem to indicate itin- 
erant moments with significant three-dimensional (3D) 
character to the magnetic ordering^i^ — far from the two- 
dimensional (2D), localized moments of the S = 1/2, 
antiferromagnetic Heisenberg behavior in the Mott insu- 
lating ground state of the cuprate parent systems^. As 
has been pointed out recently, this does not preclude 
the importance of correlation effects in the iron pnic- 
tides whose electronic ground state may simply exist in 
close proximity to the Mott phase^i^. Whether the in- 
fluence of electron-electron correlation effects within the 
underlying physics of the iron pnictides is as important 
as the known correlation effects in the cuprates stands 
as a central unresolved issue in the study of these new 
superconductors. 

Undoped, the parent compounds of the iron pnictides 
exhibit long range, three dimensional antiferromagnetic 



(3D AF) order with Fe^"*" moments aligned within the 
basal plane^'^. The onset of AF order is preceded by 
(eg. LaFeAsO)^ or accompanied by (eg. BaFe2As2)^ 
a structural phase transition from the high-temperature 
tetragonal structure to a low-temperature orthorhombic 
phase. It has been argued that the structural phase tran- 
sition is driven by the formation of long range magnetic 
orderiiSiiii, and recent work has suggested a mechanism 
of fluctuating domain boundaries between energetically 
similar magnetic states to explain this effect even in sys- 
tems where Tg > T^^. We will present an alternate 
model in this paper. As magnetism is thought to play 
a role within the formation of the unconventional, su- 
perconducting ground state in these systems^ii^ii^ii^, a 
comprehensive understanding of magnetic as well as the 
structural behavior in the iron pnictides, particularly in 
the undoped parent systems, is of considerable impor- 
tance. 

The ThCr2Si2-type bilayer AFe2As2 (A=Ba, Ca, Sr, 
...) systems, or 122-compounds, constitute one of the 
most promising avenues of experimental investigation of 
magnetism within the iron arsenides. This is due to their 
relatively large Fe^^ ordered-moments, and, most impor- 
tantly, due to the fact that it appears to be straight- 
forward to grow large single crystals of these materi- 
als. All 122 variants studied to date exhibit a struc- 
tural phase transformation from an lA/mmm tetragonal 
to Fmmm orthorhombic phase concomitant to the es- 
tablishing of an AF phase with ordered moments rang- 
ing between 0.8/i_B to 0.94/Zfliiii^ii^i2£ per Fe-site — values 
rather smaller than those predicted from band-structure 
calculation&i^i^i2^. Despite comparable moment sizes, 
there exists a significant range of T/v's among these sys- 
tems with Tn = 136K in BaFe2As2^, Tn = 205K in 
CaFc2As2^^, and T/v = 220K in SrFe2As2^^. This varia- 
tion in T/v's presumably arises from variations in the J„„ 
exchange coupling between the different compounds. De- 
tailed studies of the magnetism within the parent phases 
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of these 122-compounds, however, have yielded mixed 
results regarding the nature of both their structural and 
magnetic phase transitions. 

Neutron studies of the CaFe2As2 compound reported 
a first order magnetic phase transition at T = 173.5K 
that was coupled to a simultaneous onset of the struc- 
tural phase transition from tetragonal to orthorhombic 
symmetry Similar behavior was reported in SrFe2As2 
at T = 205K where a first order magnetic phase tran- 
sition coupled to the structural order parameter was 
observedi^i2^. This claim was supported by heat capac- 
ity measurements which implied a discontinuous tran- 
sition at the same temperatures^!, although separate 
heat capacity measurements of isostructural EuFe2As2 
reported no hysteresis in the magnetic/structural phase 
transitions and speculated the possible presence of an 
undetectably small latent heat^^. However, a combined 
Mossbauer and x-ray diffraction study claimed to ob- 
serve a continuous, second order structural phase tran- 
sition in the same SrFe2As2 and EuFe2As2 systems^^. 
Neutron studies of BaFe2As2 polycrystalline samples ini- 
tially reported a continuous onset of magnetic order be- 
low T — 143K with a small hysteresis present in the 
accompanying structural phase transition^. This was fol- 
lowed by separate single crystal studies claiming a second 
order magnetic phase transition, albeit with substantially 
renormalized magnetic properties due to Sn incorpora- 
tion within their crystal^!.. Work by a different group on 
cleaner BaFe2As2 crystals failed to observe any hysteresis 
in the magnetic phase transition's, in direct contradiction 
to separate NMR results showing a 4K hysteresis in the 
magnetic order parameter upon warming and cooling^-. 
A recent report has even claimed a hysteresis in excess 
of 20K through Tn, supporting a claim of a strongly first 
order phase transition in the same BaFe2As2 system^^. 
Clear discrepancies amongst these various measurements 
mandate continued careful studies of the phase behavior 
in these 122-systems in an attempt to resolve the under- 
lying nature of the magnetic ground state of the undoped, 
bilayer compounds. 

The fact that the magnetic phase transition in the 
BaFe2As2 system may be second order makes this com- 
pound an illuminative system for isolating the intrinsic 
nature of the phase behavior in this class of compounds. 
Additionally, the substantially reduced Tjv in BaFe2As2 
relative to its isostructural variants also makes it an in- 
teresting system for further study. The reduced T/v is 
potentially indicative of a stronger anisotropy in the spin 
behavior where 2D fluctuations within the FeAs layers 
renormalize the magnetic behavior. For these reasons, 
we chose to study the spin and structural ordering within 
the BaFe2As2 compound via neutron diffraction. 

The remainder of this paper is organized as follows: 
Section II discusses technical details of our sample prepa- 
ration, neutron scattering experiments, and bulk trans- 
port measurements. In Section III we present a prelimi- 
nary characterization of our BaFe2As2 single crystals for 
relative crystal quality comparison with existing reports 



in the literature. Section IV details our neutron diffrac- 
tion measurements probing the magnetic and structural 
phase behaviors along with preliminary measurements of 
the critical scattering above in this system. Section 
V presents our analysis of the critical behavior in this 
system and possible effects leading to discrepancies be- 
tween various studies. Finally, in Section VI we give our 
final conclusions and summarize our work. 



II. EXPERIMENTAL METHODS 

For our experiments, we grew high quality single crys- 
tals of BaFe2 As2 ( in orthorhombic notation Fmmm; a = 
b = 5.583A,c = 12.9471 at 250K) using a vertical 
Bridgman-type growth technique similar to that em- 
ployed by Morinaga et al.'^". A 657 mg single crystal was 
obtained in a sealed quartz ampoule, 12 mm in diame- 
ter, in the presence of a FeAs flux. The complete details 
of the sample growth will be published elsewhere^!. A 
366mg single crystal, the larger half of a cleaved 657mg 
single crystal, was chosen for neutron studies. Neutron 
experiments were carried out on the N5 spectrometer at 
the Canadian Neutron Beam Centre, National Research 
Council, Chalk River Laboratories. The N5 spectrometer 
was used in a triple-axis configuration with a coUimation 
of 30'— 60— sample— 52'— 72', unless otherwise indicated, 
and a fixed final energy of £'/=3.52 THz was used. The 
(0, 0, 2) reflections of pyrolytic graphite (PC) crystals 
were used for both the monochromator and the analyzer 
stages. A pyrolytic graphite filter was used to reduce 
second order contamination from the incident beam. 

The BaFe2As2 crystal was mounted in the (H, 0, 
L)scattering plane (orthorhombic unit cell notation) 
within a closed-cycle refrigerator containing He exchange 
gas to ensure efficient heat exchange between the sample 
and the cryostat. Neutron measurements were conducted 
between 3 < T < 250K in order to characterize the phase 
behavior. In this work, positions in reciprocal space will 
be referenced using reduced lattice unit (r.l.u.) notation 
in which Q{h, k,l)[A-^] = (27r/a x H,2Tr/b x K,2tt/c x 
L)[r.l.u.]. Bulk magnetization and resistivity measure- 
ments were carried out using a Quantum Design physical 
property measurement system. Our DC-magnetization 
measurements were performed on the smaller 291mg half 
of the cleaved crystal used in our neutron scattering stud- 
ies while resistivity measurements were carried out on a 
small, « lOmg, crystalline bar with dimensions 7x1.5x0.1 
mm cut from a boule from a separate growth run. 



III. SAMPLE CHARACTERIZATION 

Before presenting our neutron measurements, it is im- 
portant to discuss first the quality of the BaFe2As2 sam- 
ples studied in our experiments. Magnetization measure- 
ments of BaFe2 As2 often present the best gauge of sample 
quality; in the dc-susceptibility there emerge two salient 
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FIG. 1: (a) Bulk susceptibility measurements showing the 
DC-magnetization measured with the magnetic field applied 
parallel to the afc-plane. The dashed line denotes the AF 
transition temperature at T = 136K. (b) Neutron charac- 
terization showing a rocking curve illustrating the mosaic of 
our crystal. The resolution limited peak full width at half 
maximum (FWHM) of 0.32° demonstrates the high quality 
of our single crystals. The measurement was taken with col- 
limation of 30'— 12— sample— 12' — 72'. (c) Resistivity mea- 
surements on a small crystal of BaFe2As2 showing a linear 
decrease with decreasing temperature and a sharp drop at 
T — 135K. The drop in resistivity is indicative of the si- 
multaneous magnetic and structural phase transitions in this 
system, (d) The derivative ^ upon both warming and cool- 
ing showing a sharp peak indicative of the slope change at 
the transition temperature. The peaks during both sample 
warming and cooling overlap perfectly showing no evidence 
for hysteresis in the simultaneous magnetic and structural 
phase transition of our samples. 



features: a Curie-type upturn in the susceptibility at low 
temperatures (typically below 30K), and a sharp decrease 
in the susceptibility curve at the antiferromagnetic order- 
ing temperature. A large susceptibility upturn in the low 
temperature Curie tail is reflective of the incorporation of 
magnetic impurities, often flux, into the sample volume. 
The inclusion of impurities into the sample also broadens 
and reduces the magnitude of the discontinuity in the sus- 
ceptibility at T/y. Our dc-susceptibility measurements of 
bulk samples grown via the Bridgman technique confirm 
this qualitative statement where progressive cuts starting 
at the larger diameter end of bullet shaped boules down- 
ward to the phase-pure single crystalline tip were taken. 
These progressive cuts revealed a continual damping of 
the Curie- Weiss upturn at low temperature along with 
a sharpening of the inflection point near the magnetic 
phase transition as the single crystalline tip of the boule 



was neared. Early work on BaFe2As2 crystals grown in 
Sn-flux, which is now known to incorporate Sn impurities 
into the BaFe2As2 matrij*^^, provides a dramatic demon- 
stration of this effect where the dominant feature of the 
susceptibility is a continual Curie- Weiss-type upturn that 
persists until the lowest temperatures measured. In typi- 
cal magnetization measurements of these samples, any in- 
dication of magnetic ordering is typically obscured with- 
out subtraction of low-field data^^. 

The dc-susceptibility of our BaFe2 As2 crystal is shown 
in Fig. 1(a) where the susceptibility was measured with 
the field applied parallel to the a6-plane. A linear de- 
crease in the susceptibility from room temperature was 
observed until the AF ordering temperature of w 137K. 
Below 137K, there appears a sharp drop in the suscepti- 
bility followed by a minimal upturn in the susceptibility 
at temperatures below w 25K. The sharp phase transi- 
tion evident in our dc-susceptibility measurements com- 
bined with the excellent agreement between the magneti- 
zation curves of our bulk crystal and those of the smaller 
self-flux grown crystals in Ref. 34 demonstrates that our 
sample is of high quality and that the phase-pure section 
of the Bridgman boule had been reached. In Section IV 
of this paper, we will report from our neutron measure- 
ments that the onset of magnetic order in this BaFe2As2 
sample occurs at 136K, and this IK difference in the 
magnetic onset temperatures is likely due to differences 
in the thermometry between the neutron diffractometer 
and magnetometer sample environments. 

In order to characterize further our BaFe2As2 crystals, 
we also performed resistivity measurements in the tem- 
perature range from T = 2K to T = 300K. The results of 
our measurements are shown in Fig. 1(c). The resistiv- 
ity, p(T), continuously decreases upon cooling from 276K 
with a sharp kink in the resistivity curve observed at the 
AF phase transition. The residual resistivity ratio (RRR) 
gives p{300K) / p{5K) — 3.1, a value consistent with ear- 
lier reports'^^. The downturn in resistivity within this 
smaller sample occurs at approximately T = 135K. In 
Fig. 1 (d) the first derivative of resistivity with respect to 
temperature, is shown, and it exhibits a sharp peak 
at T = 134. 25K as expected since, in a metal, ^ scales 
like the heat capacity. The temperature difference be- 
tween the peak in the derivative of the resistivity and the 
downturn in the magnetic susceptibility is due to small 
variations in T/y between different samples where subtle 
variations in the chemical gradient in different specimens 
can alter their antiferromagnetic ordering temperatures. 
In order to perform preliminary tests for any hystere- 
sis in the magnetic and structural phase transitions in 
this system, we measured the resistivity upon both slow 
warming and slow cooling the sample through Tpf. The 
results shown in Fig. 1 (d) show no difference in the onset 
of the sharp slope-change in resistivity data upon warm- 
ing and cooling with both peaks in Sp/ST coincident at 
T = 134. 25K. This supports the notion that the magnetic 
phase transition in this material is, in fact, continuous. 

A further measure of our crystal's quality is provided 
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FIG. 2: (a) Radial scans through magnetic peaks at the (1, 0, 
L); L=l,3,5,7 positions along with the (3, 0, 3) position, (b) 
Rocking curves through the corresponding magnetic peaks in 
panel (a) . Insets in each panel qualitatively show the direction 
that the resolution ellipsoid intersects the magnetic Bragg 
position. 

by the mosaic of the crystal itself. A rocking scan through 
the (0, 0, 4) reflection at 3K, shown in Fig. 1(b), ex- 
hibits a sharp peak with full width at half maximum 
(FWHM) of 0.32°; a value determined by the resolu- 
tion of the spectrometer. This implies that the crystal 
quality is extremely high and that the crystal possesses 
a relatively small amount of strain-inducing defects and 
misalignments which may couple to and renormalize the 
intrinsic magnetic behavior. This also allows for a more 
precise determination of the critical behavior in this sys- 
tem, where the absence of multiple grains within the mo- 
saic avoids probing different crystallites with slightly dif- 
fering magnetic properties. 



IV. MAGNETIC AND STRUCTURAL PHASE 
BEHAVIORS 

A. Magnetic order parameter 

Turning now to our neutron diffraction results. Figure 
2 shows scans at T ~ 3K through the magnetic posi- 
tions Q = (l,0,L);i = 1,3,5,7 and Q = (3,0,3). Panel 
(a) of Fig. 2 displays radial scans through these mag- 
netic peaks, and panel (b) shows rocking curves through 
the same positions. Insets in each panel qualitatively 
illustrate the direction that the resolution ellipsoid inter- 
sects the magnetic Bragg peaks within the corresponding 
scans. The measured intensities are consistent with the 
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FIG. 3: Q-mapping through mesh scans about the (3, 0, 3) 
magnetic Bragg reflection with (a) the PG filter out of the 
beam and (b) the PG filter in the beam. The resulting nuclear 
scattering from the (6, 0, 6) position in panel (a) shows the 
structural twinning below Ts which is absent in the purely 
magnetic scattering from the (3, 0, 3) reflection in panel (b). 
This indicates that the moments are aligned along the longer 
a-axis in agreement with previous reports^ii^ 



published spin structure^'^^ for the 122 variants where 
the Fe-moments align within the basal plane. The mag- 
netic correlation length within the a6-plane is resolution 
limited with a minimum correlation length of 385A and 
a minimum out of plane correlation length of 350A de- 
termined with a coUimation of 30'— 12— sample— 12^— 72". 
In determining the direction that the moments lie within 
the a5-plane, the structural twinning which occurs below 
the tetragonal to orthorhombic phase transition was uti- 
lized as in Ref. 18. Maps in Q-space were performed 
with both the PG filter in the beam and with the PG fil- 
ter removed at the (3, 0, 3) position. Figure 3 shows the 
results of these scans with the PG filter out (Fig. 3(a)) 
and the PG filter in (Fig. 3(b)). In Fig. 3(a), the second 
order contamination of the beam dominates the scatter- 
ing signal through AEi neutrons scattering from the (6, 0, 
6) nuclear Bragg reflection. The distended peak shape re- 
flects the equal population of the twin (6, 0, 6) and (0, 6, 
6) domains consistent with the expected twinning upon 
cooling into the orthorhombic phase. Upon replacing the 
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PG filter in Fig. 3(b), the second order contamination 
is strongly suppressed, and the scattering that remains 
originates from the (3, 0, 3) magnetic reflection. Scatter- 
ing from both twins no longer persists, and only one peak 
appears from the single twin aligned along the longer a- 
axis. This verifies the Qj(^=(l, 0, 1) ordering vector for 
the spin structure in this system, consistent with previous 
reports^^'^^. Normalizing the (1, 0, 3) magnetic intensity 
to the (0, 0, 4) nuclear reflection yields an ordered mo- 
ment of 0.93 ± 0.06/iB per Fe^+ ion. In obtaining this 
normalization, the assumption was made that the struc- 
tural twin domain populations are equal (an assumption 
supported by our subsequent analysis of the structural 
phase transition). 

Having established that the low temperature magnetic 
structure is consistent with previous reports, we then 
examined the temperature dependence of the magnetic 
scattering at the strongest magnetic reflection at Q =(1, 
0, 3). Figure 4 shows the temperature dependence of the 
magnetic Bragg intensity upon both slowly warming and 
slowly cooling the sample. For both cooling and warm- 
ing scans, a 20 minute equilibration was allowed between 
each temperature ramp of IK with the cooling rate set 
to be 2K/min for the ramp itself. The magnetic intensity 
at the peak (1, 0, 3) position was continuously measured 
during the equilibration time resulting in one measure- 
ment per « 50mK. This was followed by a full radial 
scan for each temperature step. The temperature regis- 
tered at the sample position was coupled to the crystal 
via helium exchange gas within the sample chamber of 
the displex rendering instantaneous measurements of the 
sample temperature accurate to within lOmK. In panels 
(a) and (b) within Fig. 4, the peak intensity measured at 
the (1, 0, 3) position is over plotted with the integrated 
area of radial scans through the (1, 0, 3) position as 
the sample is respectively cooled and heated through the 
magnetic phase transition. The excellent agreement be- 
tween the temperature evolution of the peak value mea- 
sured at Q =(1, 0, 3) with the integrated radial peaks 
at discreet temperatures demonstrates that the measure- 
ments of the magnetic ordering at the peak position are 
not effected by any extrinsic effects such as temperature 
inhomogeneity within the sample and that there is no 
effect from a broadening or shift in peak position upon 
cooling/ warming through the phase transition. The tem- 
perature dependence of the Q =(1, 0, 3) peak intensity 
can therefore be used as a valid measure of the magnetic 
order parameter squared. This allows for a higher tem- 
perature resolution near the phase transition. 

We will consider first the integrated intensities of the 
(1, 0, 3) magnetic peak upon cooling and warming as 
shown in Fig. 4(c). There is no resolvable hysteresis 
and the onset of magnetic order is continuous. This is 
illustrated clearly in Fig. 4(d) where the peak intensities 
upon cooling and warming are over plotted. Fitting the 
temperature dependence of these peak intensities with a 
simple power law <j){TY oc (1 — gives an AF or- 

dering temperature of Tm — 135.95 ±0.05K upon cooling 



and Tjv = 136.02 ± 0.09K upon warming. As expected, 
they agree within the errors. Given that there is no evi- 
dence for hysteresis in the magnetic ordering in our data, 
we then fit the combined warming and cooling data to 
same power law relation giving T/y — 135.98 ±0.04K and 
13 = 0.103 ± 0.018. This is plotted as a dotted line in 
Fig. 4(d). This (3 value is close to the value, 1/8 ex- 
pected for the universality class of the two dimensional 
Ising model^ and is consistent with earlier reports by 
Matan et ali^. Upon narrowing the range of tempera- 
tures fit above and below the transition, the fit (3 value 
systematically increases toward the 2D Ising value. Re- 
analyzing (f){TY with a fixed /? = 1/8 renders the fit 
displayed as a solid line in Fig. 4(d) with a refined 
Tjv — 136.23 ± 0.025K, a value within the error of the 
result from our power law fit using a free (3 parameter. 
Given the close statistical proximity of the convergence 
of the fit (3 value to the 2D Ising exponent in the criti- 
cal temperatures closest to the phase transition and the 
fact that the constrained 2D Ising fit diverges from the 
results of the refined (3 = 0.103 fit only far from the 
transition temperature, it appears that the magnetic or- 
dering is approximated reasonably well by the 2D Ising 
exponent. It is also worth noting that the large range in 
reduced temperature over which the power law fit mod- 
els the development of the magnetic order in Fig. 4(d) 
is often a feature of other low dimensional antiferromag- 
netic systems such as K2NiF42& (see Appendix for further 
discussion) . 

To test the validity of the power law fit, we performed 
further fits probing a range of temperatures about T^y 
and a range of reduced temperatures, i, below the phase 
transition. First, we observed that, upon screening the 
temperatures closest to the phase transition and refitting 
the data, the extracted critical exponents were smaller 
than the /3 — 0.103 rendered from the full fit of the 
data. For example, fits considering only T < 130K yield 
a 2/3 = 0.171 ± 0.012 while fits considering T < 134K 
generate a 2/3 = 0.194 ± 0.006. This implies that the 
smallness of the critical exponent is intrinsic and not due 
to a sharpening induced by an unseen first order jump 
in the order parameter close to Tjv. Additionally, we 
repeated the power law fits in a shortened range of re- 
duced temperatures about the transition, and the results 
of these fits are plotted in Fig. 5 as solid yellow lines 
through the peak intensities of the (1, 0, 3) magnetic 
peak as a function of temperature. Fits within t ~ 0.10 
resulted in a Tat = 136.26±0.18K and 2/3 = 0.224±0.015, 
while fits within t « 0.05 result in a = 136.16 ±0.21K 
and 2/3 = 0.212 ± 0.021. For both fits, only data below 
T = 135. 5K were included. These (3 and T/v values agree 
with the values rendered from the full fit to the data 
within the statistical errors and therefore demonstrate 
that our fit results are robust. 

In order to constrain the range over which we can pre- 
clude the existence of a weakly first order phase tran- 
sition in this system, we performed a linear fit to the 
magnetic order parameter squared at the temperatures 
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FIG. 4: Measurements of the magnetic order parameter upon 
warming and cooling at the (1, 0, 3) reflection. Both the peak 
intensity measured at Q =(1, 0, 3) and the radially integrated 
(1, 0, 3) peak area are plotted upon both (a) cooling and (b) 
warming through Tjv. The perfect overlap between the peak 
intensity and integrated area measurements indicates that the 
peak intensity can be safely used as a measure of the mag- 
netic order parameter, t/)^. (c) Integrated (1, 0, 3) peak areas 
overplotted as a function of cooling and warming. There is no 
sign of hysteresis and the onset of magnetic order is seemingly 
continuous, (d) Peak intensities at Q =(1, 0, 3) overplotted 
upon both cooling and warming the sample. The dashed line 
represents the results of the power law fit described in the 
text with P — 0.103 while the solid line displays the results of 
a power law fit with l3 fixed at 1/8. 

closest to the phase transition. This is shown within the 
inset of the top panel in Fig. 5. The break in the data 
between T — 135. 8K and 136. IK simply constitutes a 
small temperature range where equilibrated data were 
not collected and not a discontinuity in the order param- 
eter. The linear fit in this region represents the linear 
trade-off which must occur between phases during the 
onset of a first order phase transition, and, from the fit 
plotted within the inset of Fig. 5, our data allow us to 
preclude such a linear trade-off below 135. 5K. Therefore, 
any first order nature to the phase transition is intrin- 
sically weak and buried within the small temperature 
interval 135.5 < T < 136. OK; however it would seem 
to have little consequence to the overall phase transition 
that appears to be fluctuation driven (as discussed fur- 
ther in Section IV C). 



B. Structural phase transition 

Simultaneous to the onset of long range magnetic or- 
der, we observe a splitting of the (2, 0, 0) structural peak 
at Tg = 136. OK indicative of the structural phase transi- 
tion from 14 1 mmm to Fmmm symmetry and consistent 
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FIG. 5: Continued analysis of the magnetic order parameter 
in different regimes of reduced temperature. Top panel shows 
a power law fit to the Q =(1, 0, 3) peak intensity data in the 
reduced temperature regime t ~ 0.10 while the bottom panel 
shows a power law fit to the data in the t « 0.05 regime. For 
both fits, only data below T — 135. 5K were included. Power 
law fits are denoted by the solid yellow lines in both panels. 
The inset in the upper panel replots the magnetic phase tran- 
sition zoomed in to within ±1K of the phase transition. A 
linear fit to the data between 135. 6K and 136. 3K is plotted 
as a solid line through the data and illustrates the maximum 
possible range for a linear trade-off between the paramagnetic 
and ordered magnetic phases required by a first order phase 
transition. 



with other observations of simultaneous magnetic and 
structural phase transitions in 122 variants^ii^ii^. Due to 
resolution constraints and the intrinsically small struc- 
tural distortion, the splitting of the (2, 0, 0) reflection 
could not be resolved into two distinct peaks; however, 
a clear broadening could be observed. The broadened 
line shape may be effectively fit by assuming two resolu- 
tion limited Gaussian peaks split symmetrically from the 
Q =(2, 0, 0) position. This reflects the simply convolved 
projection of the resolution function along the a-axis with 
two intrinsically small nuclear peak widths approximated 
by delta functions. The quality of the fits below and 
above Ts is shown through the solid lines in Fig. 6(a), 
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where the broadened peak at T = 98K is fit by two Gaus- 
sian hne shapes spUt by a parameter 5 at (2±(5, 0, 0) and 
at r = 140K where the spHtting is effectively i5 = . The 
parameter 5{T) is directly representative of the structural 
order parameter via the relation | — . The structural 
order parameter squared (^(T)^) is then plotted in Fig. 
6(b) for both slow cooling and warming of the sample, 
using cooling/warming rates identical to those previously 
described in the measurement of the magnetic order pa- 
rameter. It is immediately apparent that there exists no 
observable hysteresis within the structural phase transi- 
tion with both slow cooling and warming scans yielding 
an identical Tg — 136. OK concomitant with the onset of 
magnetic order. 

The smooth, non-hysteretic evolution of the structural 
order parameter through Tg in Fig. 6(b) suggests a con- 
tinuous transition coinciding with the onset of long-range 
antifcrromagnetic order within this system and consistent 
with the theoretical understanding of a structural distor- 
tion coupled to the magnetic ordering2iiiii^. To explore 
this further, we compared the temperature evolution of 
both 0(T)^ (the magnetic order parameter squared) and 
5{TY upon both warming and cooling as shown in Figs. 
6(c) and (d) respectively. Both panels (c) and (d) dis- 
play a remarkable agreement between the square of the 
magnetic and structural order parameters suggesting a 
strong biquadratic coupling term e^0^ between the two 
order parameters. This contrasts with the behavior ex- 
pected of the structural distortion as a secondary order 
parameter driven by the magnetic order which functions 
as the primary order parameter. In that case there is 
an expected linear-quadratic coupling of the strain and 
staggered magnetization (cx £0^). This makes the unex- 
pected agreement between the <5(T)^ and 4>{T)'^ parame- 
ters particularly interesting, and the implications of this 
coupling will be discussed in Section IV. 



C. Critical scattering 

Given the apparently continuous nature of the mag- 
netic order parameter within Figs. 4 and 6 and the lack 
of any observable hysteresis in the AF ordering tempera- 
ture, we next searched for the presence of magnetic crit- 
ical scattering near Tjv. The analyzer was left in for 
these measurements due to a gain in the signal to back- 
ground ratio. Hence any critical scattering observed does 
not represent an energy integrated response, but rather 
only the quasi-elastic critical scattering integrated over 
the spectrometer's energy resolution {AE = 1.14 meV 
FWHM). The results of our critical scattering measure- 
ments are shown in Fig. 7. The scattering about the 
strongest magnetic reflection, Q =(1, 0, 3), was mea- 
sured in the temperature range immediately above Tn 
revealing a weak critical scattering signal extending up 
to T « 150K (shown in the shaded area of Fig. 7(a)). 
Each point in Fig. 7(a) represents the integrated area 
of radial scans through the (1, 0, 3) position, and select 
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FIG. 6: The evolution of the structural order parameter as 
a function of temperature, (a) Longitudinal H-scans through 
the (2, 0, 0) nuclear Bragg reflection both below and above 
the structural phase transition Ts. The solid lines are the 
results of fits of two Gaussians to the broadened peak below 
Ts with the splitting, 5{T), determined as described in the 
text, (b) (5(r)2 plotted upon both warming and cooling. The 
onset of the structural splitting shows no hysteresis, identical 
to the development of the magnetic order in this material. 
(5(r)^ is overplotted with the integrated area of the (1, 0, 
3) magnetic reflection as a function of (c) warming and (d) 
cooling showing a remarkable tracking between the structural 
order parameter squared and the magnetic order parameter 
squared. The exact tracking within the errors verifles previous 
suggestions of a strong coupling between the onset of magnetic 
order and the structural lattice distortion in this material^ 



scans at temperatures above T/v are plotted in Fig. 7(b). 
The residual peak above T = 150K results from remnant 
second order contamination in the incident neutron beam 
scattering from the nuclear (2, 0, 6) reflection. This was 
determined by the lack of temperature dependence of the 
remnant (1, 0, 3) signal upon warming from T = 150K 
to T = 250K. 

The dimensionality of the critical signal was measured 
via scans in Q-space through the (1, 0, 3) position prob- 
ing spin correlations along the a-axis and correlations out 
of the afc-plane, along the c-axis. The results of these Q- 
scans along H and L are plotted in Figs. 7(c) and (d). 
In order to isolate the critical component of the peaks, 
the residual A/ 2 nuclear peak was subtracted off using 
data collected at T = 148K as the nuclear background. 
Sharp, resolution-limited peaks centered at the Q =(1, 
0, 3) position in both H and K scans above T/v reveal the 
critical scattering to be fully three dimensional. However, 
the widths of peaks in both H and L— scans were reso- 
lution limited at all temperatures precluding any anal- 
ysis of the magnetic correlation length as a function of 
temperature. Additionally, scans at non-zero offset 
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from L = 3 revealed no evidence for a substantial two 
dimensional component to the critical scattering. This, 
however, does not preclude the existence of such two di- 
mensional scattering where the signal may be too weak 
to be resolved due to the relatively small sample volume 
of our crystal. Proper two-axis measurements integrat- 
ing over a large range of energies on a larger sample or 
sample array would provide a more conclusive probe of 
any putative two dimensional critical scattering in this 
system. 

The existence of critical scattering above the phase 
transition supports the observation that the phase tran- 
sition is second order; however, the strongly three dimen- 
sional nature of the critical scattering above T/v contra- 
dicts the behavior expected for phase transitions within 
the 2D Ising regime. Then, at the minimum, there must 
be a crossover from 2D to 3D critical behavior near Tjv. 
The three dimensional nature of the critical scattering 
can be reconciled with the anisotropic spin wave disper- 
sion reported by Matan et al.^" with a c-axis spin wave 
velocity of Vc = 57meVA relative to the in-plane spin 
wave velocity of Vab = 280meVA. The two dimensional 
inelastic scattering above T/v observed by Matan et al;^ 
does not contradict our observation of three dimensional 
critical fluctuations driving the magnetic phase transi- 
tion. Indeed the observed behavior is similar to that of 
FeClz^, where the critical behavior was observed to be 
three dimensional despite a factor of 20 difference be- 
tween in-plane and out-of-plane exchange couplings. 



V. DISCUSSION 

The first issue to be addressed is the direct discrepancy 
between our neutron measurements revealing a seemingly 
second order magnetic phase transition in BaFe2As2 and 
those recently reported by Kofu et al.^— which claimed 
to observe a strong first order magnetic behavior. One 
key factor which can potentially alter the magnetic phase 
transition in the iron pnictides is the coupling of the mag- 
netic order to elastic strain within the system. It has long 
been known that through coupling to elastic strain or 
lattice deformation a second order magnetic phase tran- 
sition can be driven first order—. Until recently, this 
effect has not been taken into account in the iron pnic- 
tides; however, recent work by Barzykin and Gorkovii 
showed explicitly the ability of an external homogeneous 
deformation to alter dramatically the resulting magnetic 
phase behavior in the iron arsenides. The application of 
a lattice deforming strain can transform the order param- 
eter through modifying the degree of Fermi surface nest- 
ing responsible for the formation of the SDW state, and, 
given sufficient strain, the Landau functional is modified 
to exhibit first order behavior. However Barzykin and 
Gorkov^ also predicted that the orthorhombic transi- 
tion itself imposes an intrinsic quadratic magnetostric- 
tion term in the free energy which is alone sufficient to 
drive the transition to first order. Looking at the large 
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FIG. 7: Critical scattering above Tm at the (1, 0, 3) magnetic 
reflection, (a) Radially integrated area of the (1, 0, 3) mag- 
netic peak as a function of temperature. The critical scatter- 
ing above Tjv is denoted by the shaded area. The dashed line 
denotes the level of A/2 contamination in our measurements 
which can be treated as background. Radial scans at select 
temperatures are plotted in panel (b). The large A/2 peak ap- 
parent at 148K renders this measurement of critical scattering 
only preliminary. Q-scans through (1, 0, 3) along both the 
H and L directions in reciprocal space are plotted in panels 
(c) and (d) respectively. The nuclear background scattering 
from the 148K A/2 peak was subtracted off in order to isolate 
the magnetic portion of the peaks. Sharp, resolution-limited 
peaks in both in plane and out of plane scans demonstrate 
that the critical scattering observed is fully three-dimensional. 
The solid lines in panels (b)-(d) are the results of Gaussian 
fits to the data. 



variance of behavior between our measurements of mag- 
netic ordering in BaFe2As2 and those performed by Kofu 
et ali^, it is reasonable to infer that this intrinsic stric- 
tion induced upon cooling through Ts does not repre- 
sent the dominant strain in this system. The samples in 
both our study and that by Kofu et ali^ should possess 
the same intrinsic magnetostriction terms upon cooling 
through the phase transition yet they yield vastly differ- 
ent outcomes. This leads us to speculate that external 
strain imposed while mounting the crystal or sample de- 
pendent internal strain fields formed during growth due 
to imperfections or impurities plays the dominant role in 
modifying the phase transition. 

Typical studies of 122 single crystals, consist of exper- 
iments on relatively thin, plate-like crystals with thick- 
nesses on the order of 0.1-0.2 mm and much larger sur- 
faces parallel to the a6-plane (often w 5mm x 5 mm or 
greateriiSii^) . This crystal geometry is especially suscep- 
tible to external strain upon mounting and possibly to 
the formation of internal strain fields due to dislocations 
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during crystal nucleationSS,. One notable effect of strain 
in these smaller samples is often manifested in the large 
asymmetry in their twin structural domain populations 
below Ts, where, by symmetry, there is no energetic dif- 
ference between domains aligned along the orthorhombic 
a- and fe-axes^^. 

Another factor which must be considered is the role 
of potential impurities within our BaFe2As2 sample. Im- 
purities and imperfections have been observed to renor- 
malize a first order phase transition into a continuous 
phase onset"^^ . One hallmark of this effect is a broadened 
critical regime remnant of the smeared intrinsic first or- 
der transition. Figures 3 and 4, however, show a sharp 
transition in our crystal which reaches half of its satu- 
ration value below Tjv at T'50% < IK. Comparing this 
to previous studies of magnetic order in BaFe2As2 re- 
veals that in fact our transition is much sharper than 
previous studies reporting discontinuous behavior where 
the T^a% ~< 20K^ in powder samples and T'50% ~< 7K 
in single crystals^^. This fact, combined with the sharp 
mosaic and clean transport properties of our sample, pro- 
vides compelling evidence that the continuous transition 
we observe is not simply an impurity-induced renormal- 
ization of an otherwise first-order transition. We there- 
fore argue that strain or other extrinsic effects play the 
dominant role in renormalizing the phase transitions re- 
ported in previous studies of BaFe2As2 where the mag- 
netic phase transition is broadened relative to the phase 
transitions observed in our work. It is, however, not pos- 
sible to preclude the existence of a weak first order tran- 
sition within 0.5K of Tjv. The presence of hysteresis in 
such a transition is difficult to detect given the small (3 
value of the power law behavior; however existing reports 
of hysteresis in excess of 4 K^^i^^ would have been easily 
detected in our measurements. 

Additionally, a small first order hysteresis can be sub- 
tly smeared over a narrow temperature range in the crit- 
ical regime generating a seemingly continuous response 
that is difficult to discriminate against. As discussed in 
Section IV A, we performed additional fits in order to 
ensure that the small critical (3 extracted from our mea- 
surements is not artificially low due to the influence of 
such a smeared first order transition through screening 
those temperatures closest to the magnetic phase onset. 
In all cases, as the range of temperatures fit moved fur- 
ther away from T^r, the subsequently refined (3 became 
progressively smaller than the /3 = 0.103 extracted from 
the full fit of the data. This strongly suggests that our 
observation of a sharp critical exponent consistent with 
the (3 — 1/8 of the 2D Ising model is not an artifact from 
an unseen first order phase transition artificially sharp- 
ening the onset of the magnetic phase. 

A key result of our study is the unexpected behavior 
of the structural phase transition in BaFe2As2. The si- 
multaneous onset of both the structural and magnetic 
phase transitions in this system is consistent with the 
current experimental pictures of structural and magnetic 
order in the undoped 122 compounds; however, the in- 



creasingly prevalent speculation that the magnetic order 
drives the structural phase transition in BaFe2As2 is con- 
tradicted by our analysis of the relationship between the 
structural and magnetic order parameters. Within the 
scenario of a secondary structural order parameter driven 
by the magnetic order through magnetoelastic coupling, 
the Landau free energy expansion for an order-disorder 
phase transition in an Ising system is of the form'^S: 
Fme = Di(j)'^e + O.SGi^^e^ -t- ... The key term in this 
expansion is the linear-quadratic coupling term of the 
structural and magnetic order parameters respectively. 
Given this term, the square root of the structural order 
parameter is expected to track the magnetic order param- 
eter as a function of temperature. Instead, as shown in 
Figs. 6(c) and (d), the structural order parameter tracks 
the magnetic order parameter directly. This strongly sug- 
gests a biquadratic coupling term between the structural 
and magnetic order parameters in the free energy of this 
system and that the structural order parameter enters 
the free energy expansion as a primary order parameter, 
a{T). 

The simplest free energy expansion of a system with 
coupled magnetic and structural order parameters is: 
Ftetra = ^\T) + a{Tf + 2d<j,{Tfa{Tf + a\T) + 
(j)'^{T) + ..M'^"^ where </'(T) is the magnetic order pa- 
rameter, a(T) is the structural order parameter, and d 
is the coupling parameter between the two. This free 
energy functional implies the existence of a tetracritical 
point in the phase diagram of the system. Tetracritical 
points have been previously observed within the phase 
diagrams of the high-Tc cuprates through a biquadratic 
coupling of the spin density wave and the superconduct- 
ing order parameters^. Theoretical treatment of this 
coupling generated exotic phase diagrams as a function 
of both doping and field tuning^i^. In the proximity of 
the tetracritical point, a mixed phase of coexistence be- 
tween the two coupled phases a{T) and (^(T) exists below 
the intersection of their two second order phase transition 
lines, and the extent of this mixed phase within the phase 
diagram is a function of the coupling strength d. The 
simultaneous onset of both the magnetic and structural 
phase transitions within BaFe2As2 strongly suggests that 
the system lies very close to the tetracritical point within 
its phase diagram. While doping potentially changes 
many parameters simultaneously, the onset of the cou- 
pled phases should subsequently split as the system is 
tuned away from the tetracritical point regardless of the 
microscopic details of the system. In fact, such a split- 
ting has been recently reported in the phase diagram of 
Ba(Fei_a;Co2,)As2'*^ which further supports the picture 
of a tetracritical point in this system. Furthermore, this 
picture of two coupled primary order parameters in the 
close proximity of a tetracritical point supports the no- 
tion that the phase transitions are second order. For an 
intrinsically first order transition in the transition region 
there must be a linear trade-off between one phase and 
another, and this could only happen in our case very close 
to Ttv, which is not what is observed. 
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The seeming contradiction between the /3 = 1/8 two- 
dimensional Ising exponent modehng the critical regime 
of BaFe2As2 and the three dimensional critical fluctua- 
tions apparently driving the phase transition can be un- 
derstood when the known asymmetry of the dimensional 
crossover within 2D Ising phase transitions is taken into 
account. For a quasi- two dimensional Ising magnet with 
a finite coupling R between the layers, there must ex- 
ist a crossover to three dimensional magnetic behavior 
at temperatures close to phase transition. The observa- 
tion of three dimensional critical fluctuations above Tn 
combined with the 2D Ising exponent modeling the de- 
velopment of magnetic order indicates such a crossover 
from two dimensional to three dimensional behavior in 
BaFe2As2. If we assume that a measurable crossover 
to three dimensional phase behavior manifests itself at 
t sa 0.050 in reduced temperature above T/y, then be- 
low T/v such a crossover would persist to only t « 0.006 
in reduced temperature^i^. This is due to the known 
asymmetry in crossover regimes of the two dimensional 
Ising magnet close to the phase transition. The small 
crossover temperature range below T/y is insufficient to 
modify the overall power law fit to the critical behav- 
ior and the three dimensional behavior near the critical 
regime is effectively washed out by the much broader 
range of temperatures at which two dimensional behav- 
ior persists. This conjecture would be confirmed by the 
observation of a crossover from 3D to 2D critical fluctu- 
ations above Tjv. 

VI. CONCLUSION 

In summary, we have studied both the magnetic and 
structural phase transitions within a large, high quality 
single crystal of BaFe2As2. Our results reveal simulta- 
neous, continuous magnetic and structural phase tran- 
sitions at Tjv = 136. OK. The lack of hysteresis in the 



order parameter onset temperatures and the detection 
of 3D critical scattering suggest that the phase transi- 
tion is second order, although it is impossible to pre- 
clude a weak first order transition whose discontinuity 
is within 0.5K of T^r. Our measurements demonstrate 
that within bulk single crystal specimens of BaFe2As2 the 
substantially reduced strain from crystal mounting and 
improved crystal quality allows for the intrinsic phase 
behavior of BaFe2As2 to be resolved. This is supported 
by consistent results between our bulk sample measure- 
ments and those by Matan et al32. whose sample was 
grown via similar methods. We also report biquadratic 
coupling between the structural and magnetic order pa- 
rameters which strongly suggests that BaFe2As2 is in 
the immediate vicinity of a tetracritical point. Further 
work exploring the detailed doping evolution of the mag- 
netic and structural phase transitions in this system is 
required to explore this scenario fully. The critical be- 
havior of the magnetic order parameter is well fit by 
the two-dimensional Ising critical exponent; however, the 
critical scattering immediately above T/v is observed to 
be purely 3D. Further experiments, with larger sample 
volumes, are required to make a definitive statement re- 
garding the presence of any 2D component to the critical 
scattering as would be implied by the 2D Ising behavior 
observed below Tpf. 
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